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Abstract. We derive a (l+l)-dimensional nonlinear evolution equation (NLE) which 
may model the propagation of high-frequency perturbations in a relaxing medium. As 
a result, this equation may possess three typical solutions depending on a dissipative 
parameter. 
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I Nonlinear dynamics may be a topic of interest in various fields of science and 
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engineering. A lot of problems arising in science and engineering may be modelled as a 
dynamical system. As an illustration, there is the Van der Pol equation [1] given by 

d 2 u ^ ON du 



_- M l-^)_ +1 ^ = 0, (1) 

found in nonlinear circuit theory [2] . The quantity u is a physical observable depending 
on the time t; fi and ujq are constants. There is also the Van der Pol-Duffing equation [3] 
given by 

d/^u du 

— - fi(l-u 2 )— + uj 2 u + uu 3 = 0, (2) 

which may model the optical bistability in a dispersive medium [4]. The additional 
quantity v is a cubic parameter. 

Besides, nonlinear phenomena may be also described by nonlinear partial 
differential (NLPD) equations such as the well-known Vakhnenko equation [5] given 
by 

Uxt + 7^{u 2 )xx + u = 0, (3) 
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arising in relaxing media as a model equation of propagation of high-frequency 
perturbations. Subscripts denote partial differentiation with respect to time t and space 
x. Equation ([3]) has been subject to many investigations ( [6-9] and references therein) in 
recent years. One typical class of solutions to NLPD equations are the so-called solitons 
arising as the result of balance between nonlinear and dispersion effects. Higher order 
solitons may be also found in higher order NLPD equations. One important question 
that may be pointed out is whether higher order solitons may survive in higher order 
Vakhnenko equation [5]. 

In the present letter, we consider a barotropic medium p = p(p, A) under relaxation. 
The quantities p and p denote pressure and mass density, respectively, and A is an 
additional parameter. Then, we derive a novel (l+l)-dimensional NLE model equation. 
We discuss the different soliton solutions to this (l+l)-dimensional NLE equation. 

Recently, Vakhnenko [5] has derived a dynamic state equation using an expansion 
of the specific volume V as power series of the small perturbations p' <C p with accuracy 
o{p' 2 ). The quantity po is the pressure related to the unperturbed state. Performing 
this expansion to accuracy o(p' 3 ), the following dynamic state equation may be found 

r {p' xx ~ Ptt/vj + a f (p' 2 ) u + a f (p' 3 ) a ) t + p' xx - p'Jv 2 + a e (p' 2 ) a + a e (p' 3 ) u = 0, (4) 

where r is the relaxation time, = ■^^ L \ p = Po , a, = g^^r|p=p anc ^ v i ma Y stand 
for velocities of the relaxation processes defined as v 2 = For high-frequency 

perturbations, that is i = f and p' = p'j = p'(p, 1) and for low-frequency perturbations, 
that is i = e and p' = p' e = p'(p, 0). Following the ref. [5], equation (TJJ may be analyzed 
by means of the multiscale method [10, 11] by introducing a small parameter e = tuj 
where uj is the frequency of the wave perturbation. 

Using a dispersion relation of the form u = v e k + j(3 e k 2 — 7 e /c 3 where j 2 = —1 for 
the linearized equation (0J, in the case of low-frequency perturbations, that is tuj <S 1, 
a (l+l)-dimensional NLE equation may be derived as follows 

Pt + v ePx + a e v 3 e (p 2 ) xx + a e v\ (p 3 ) x - (3 e p' xx + 7eP xxx = 0, (5) 

with 

a = §M-« 3. * = ^W-3W-*i). («) 

The nonlinear terms have been reconstructed in agreement with the initial equation. 
This equation (JHJ) may be viewed as a modified Korteweg-de Vries-Burgers (mKdVB) 
equation. Recently, Fu et al. [12] have studied the specific case a e = and (3 e = 
(mKdV equation) by constructing breather lattice solutions. Similar procedure may be 
discussed for equation (J5j) in order to find out other kind of breather lattice solutions. We 
may actually think that equation (jSJ) may deserve further interests from the viewpoint 
of investigation of propagation of localized and periodic waves. 

In the case of high-frequency perturbations, that is tuj ^> 1, performing the 
dispersion relation to vJ 2 uj = k 2 + jf3fk — 7/, in agreement with the initial equation, 
one may get the following (l+l)-dimensional NLE equation 

PL - vj 2 p' tt + Ufv) (p' 2 ) xx + a s v) (p' 3 ) xx + Pfp' x + lf p' = 0, (7) 
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where 

2 2 4 4 

V f — V g vj — vl 

" f = 7/ = 2r 2 ^' (8) 

standing for dissipative and dispersive parameters. In order to investigate the equation 
([7]), it seems useful to consider the following accuracy 

d 2 x -vfd^2d x (d x + vj 1 d t ). (9) 

We consider two interesting cases: atj = and aj ^ 0. 

(i) First case: ctf = 0. 

Equation ([7]) may be reduced to 

Uy V ~ ^{U 3 )yy + CXlLy ~ U = 0, (l0) 

up to the following transformations 




V= ) J%(v f t-x), V=^v f t, « = fr^, P=^=- (11 

Without dissipative a-term, equation ffTU]) may be observed as the Schcifer- 
Wayne short pulse (SWSP) equation [13] which has been subject to many recent 
investigations [14-20]. This SWSP equation may have a variant form given by 

u xt + l(u 3 ) xx + u = 0, (12) 
o 

up to the transformations x^jx, y—>jy, t—^jt and u—>-ju, j 2 = —1. Further 
interests ought to be paid to equation (TlTJT) which may have many applications 
in soliton theory and nonlinear optics. In particular, its extension to a complex- 
valued SWSP equation [19] with dissipative a-terms may be worth investigating 
alongside the effect of the dissipative parameter a on the different solutions. Indeed, 
extending u to a complex-valued quantity Q [19] in order to get the following 
equation 

Q yn - l -{\Q\ 2 Qy) y + aQ y -Q = ^ (13) 
equation ( JTBl may be transformed into the following system 

Qla ~ = {Z* + Z T )Q r - + 

QZ ~ Q™ = (Z* + Z T )Q im - a(QT + QT), (14) 

z aa - z TT = -Q r (Q r a + q T ) - Q im {QT + Q7), 

up to the following transformations 



1 



T 



y = -[X + -l QQ*dT')+n, V = T, (15) 



where \i is an arbitrary constant, T — h(a — r) and X = — |(cr + r). Looking for 
soliton solutions with the boundary conditions \Q\ — > 0, Z — > a/2 as \a\ — > oo, 
equation (fl4l may be bilinearized as in ref [19] according to the Hirota's method 
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[21,22], and soliton solutions to equation (TT4"|) may be easily derived. Thus, deriving 
the dispersion relation which may obviously be expressed in terms of the dissipative 
parameter a, it is then possible to discuss the soliton solutions of equation ( Tl3l) with 
respect to a. As a result, a soliton solution Q may be given by 

Q = Asech{$ r ) exp(^ m ), (16) 

where d = ko~ — lot + $o, being a contant parameter, and 

A = A{k r + u r ). (17) 

The physical complex-valued quantities k and ui standing for wave number and 
angular frequency, respectively, may satisfy the following dispersion equation 

k 2 - cu 2 + a(k - to) - 1 = 0, (18) 

from which k and to may be expressed in terms of the dissipative parameter a. 
Thus, discussing the soliton solutions expressed in terms of Q r and Q im vs y, one 
may expect to find loop-, cusp- and hump- shaped solitons. For some convenience, 
we do not go further with these developments. We shall focus our attention to a 
further novel equation below that may be also of great interests. 

(ii) Second case: atf ^ 0. 

Equation §Z§ may be reduced to 

*(*♦■*♦**)•♦• - 

up to the following transformations 



Without the dissipative term and d y ( ^d y uj-texm, equation ffT§j) may be reduced 
to the well-known Vakhnenko equation [5]. Without the dissipative term and 
d y (w<9yM)-term, equation (fT9l) may be reduced to (|T2|) . Performing variable 
transformations, we introduce new independent variables £ and ( as follows 

V = (+ I (u + lu 2 )d£ + y , = (21) 



2 

where y is an arbitrary constant. Then, equation fTTTJT) is reduced to 

+ + <£>w = 0, (22) 

where 

1+ / u c (1 + u) di'. (23) 

J —oo 

Defining another independent variables a and r as follows 

£=^-t), C=~(^ + r), (24) 
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equation (I22p is transformed to 

u ua — u TT = (pu — a{u a + u T ). (25) 
Moreover, using the ansatz 

<p = -Z c = (Z a + Z T ), (26) 
we get the following coupled equations 

U aa - Utt = (Z a + Z T )u - 0>(u a + U T ) , 

Z - Z TT = -u{u c + u T ) - (u a + u T ). 

This system may be closely related to that described by Kakuhata and Konno [23] 
while investigating the loop soliton solutions of string interacting with external 
field. Thus, the other physical meaning of equation (fl9l) is pointed out. This may 
be useful in constructing the soliton solutions to equation ( Fl9l) . Thus, in order to 
find a soliton solution, we consider the following boundary conditions 

it — > 0, Z — > a/2, as a — > — oo. (28) 

We may consider the following settings [23] 

u=%, Z=ha + r) + 2(d T -d a )lnF. (29) 

Equation fl271) is then bilinearized as follows 

[pl -D 2 T + a(D a + D T f - 1) (F ■ G) = 0, 

1 (30) 
(D a - D T f (F-F)-- (G 2 + 2GF) = 0, 

where D a and D T denote Hirota operators [21,22]. Expanding F and G in a suitable 
formal power series, a soliton solution to equation (|27[) is given by 

u = 4(co> + k) 2 [tanh(^) + 1] , Z = -(a +r) — 2(u + k) [tanh(6>) + 1] . (31) 

where 9 = ka — lot + 6q, 6q being an arbitrary constant. The dispersion relation is 
given by 

A(k 2 - to 2 ) + 2a(k - u) - 1 = 0, (32) 
which may lead to the following solution {k > 0) 

k = , to = kv, (33) 

«(1 -v) + v/a 2 (l - v) 2 + 4(1 - v 2 ) 

where v is the velocity of the wave satisfying the condition — 1 < v < 1. It seems 
worth noting here that from equations (1211) . fT23l) and (TSB"]) . one may find that 
y = —Z + C, C being an arbitrary constant. In order to discuss the soliton solutions 
to equation ffl9l . it is important to consider the following relation 

d a = ~ [1 - A(lo + k)ksech 2 (9)] 8 Z . (34) 

We may pay interest to the shape of the soliton u and its momentum 7r = u a + u T . 
As a result, it comes that 
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• for a = v \_J , uz may never change sign but may be infinite at some 'singular' 
point, whereas ttz may change sign once and may be infinite at the same 
singular point. Thus, u may be monotone but may have an infinite derivative 
at this particular point, and n may have a cusp-like shape (see panels [T](a) and 

HP>)); 

, Uz may change sign twice and may be infinite at two 
singular points, whereas itz may change sign three times. Thus, u may follow 



for a G 
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a multi-valued shape with two singular points at their derivatives and tt may 
have a multi-valued profile especially a loop-like shape (see panels [flc) and 

Hid)); 



finally, for a G 



oo 



Uz may never change sign and may never take 



1— V ' 

infinite values, irz may change sign once and may always be finite. Thus, u 
may have a kink-like shape, whilst 7r may have a single- valued profile especially 
a hump-like shape (see panels [D^e) andCQT)). 

We give some illustrations of the previous discussions. Thus, we may take a velocity 
v = 0.24 to plot the different profiles. The aforementioned shapes are clearly 
depicted in figure [Q, at initial time r = 0. Particularly, for the cusp-shape, the 
dissipative parameter is given by a = 0.351648275547. 

In conclusion, the studies of the novel (l+l)-dimensional NLE equation (Tl9l including 
the Vakhnenko and the variant SWSP equations (j3J) and ( fT2|) . respectively, may have 
some scientific interests both from the viewpoint of the investigation of the propagation 
of high-frequency perturbations and from the viewpoint of the existence of stable 
wave formations. Thus, applications may be found in soliton theory, geodynamics, 
hydrodynamics and nonlinear optics, just to name a few. 
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